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ABSTRACT 


Effects  of  surface  heat  transfer  on  boundary-layer  transition 
are  analyzed  in  a  three-part  study  using  the  Saffman-Wilcox 
transition  model.  In  the  first  part  of  the  study,  model  pre¬ 
dictions  are  compared  with  experimental  data  for  cooled  and 
heated  aerodynamic  boundary  layers  on  smooth  flat  surfaces  and 
for  cooled  aerodynamic  boundary  layers  near  the  stagnation 
point  of  a  roughened  blunt  body.  Consistent  with  measure¬ 
ments,  tae  model  predicts,  on  the  one  hand,  that  heating 
destabilizes  a  smooth-surface  aerodynamic  boundary  layer 
and,  on  the  other  hand,  that  cooling  destabilizes  a  rough- 
surface  aerodynamic  boundary  layer.  Differences  between 
predicted  and  measured  transition-point  locations  are  with¬ 
in  experimental  error  bounds.  Then,  incipient  transition 
conditions  are  determined  for  a  small,  heated  hydrodynamic 
body.  Again  model  predictions  agree  with  m.easurements  which 
indicate  that  relatively  small  amounts  of  surface  heating 
have  a  strong  stabilizing  effect  on  hydrodynamic  boundary 
layers.  In  the  final  part  of  the  study,  transition  location 
is  determined  for  a  large  hydrodynamic  body;  results  indi¬ 
cate  that  large  surface  heating  rates  are  not  substantially 
more  effective  than  smaller  rates.,; 
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NOTATION 


SYMBOL 


DEFINITION 


e 

F(n) 

J 

k 

^1 

P 

Pr,  , 


Pr, 


r 

R 


o 

Re^, 

Reg 

Rey 

s 

"l 

t 

T 

'^t 

00 

T» 

u,v 

u<„ 

V 


n 


Re 


Pressure  coefficient,  (p-Poo)/^  ^ 

Specific  heat  at  constant  pressure 
Specific  turbulent  energy 

?’arsons-Goodson— Golds chmled  shape  parameter 
Nondimensional  self-similar  velocity  profile 
0  for  two-dimensional  flow;  1  for  axisymmetric  flow 
Roughness  height 

Parsons-Goodson-Goldschmied  shape  parameter 
Pressure 

Laminar  and  turbulent  Prandtl  numbers 
Surface  heating  rate;  average  value  of 
Radial  distance  from  symmetry  axis 
Parson-Goodson-Goldschmied  shape  parameters 
Empirical  constant 

Reynolds  number  based  on  arclength,  plate  length 

Reynolds  number  based  on  momentum  thickness 
Incipient  transition  Reynolds  number  based  on  V 
Arclength 

Parsons-Goodson-Goldschmied  shape  parameter 
Parsons-Goodson-Goldschmied  shape  parameter 
Temperature 

Freestream  total  temperature 
Freestream  turbulence  intensity 
Velocity  component  in  x,y  direction 
Boundary-layer-edge  velocity 
Freestream  velocity 
Body  volume 
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NOTATION  (cont.) 


SYMBOL 


DEFINITION 


x,y 

m 

a, a* 

n 

e 


K 

x,x« 

A 


y 

V 


p 

a, a* 
D 


Distance  parallel  to,  normal  to  body  surface 

Parsons-Goodson-Goldschmied  shape  parameter 

Empirical  parameters 

Empirical  constants 

Similarity  variable 

Momentum  thickness 

Karman ' s  constant;  thermal  conductivity 
Empirical  parameters 

Modified  Polhausen  pressure  gradient  parameter 

Molecular  viscosity 

Kinematic  viscosity 

Fluid  density 

Empirical  constants 

Specific  dissipation  rate 


subscripts 
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Boundary  layer  edge 
Transition  point 
Surface 
Freestream 
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1.  INTRODUCTION 


In  recent  years,  substantial  progress  has  been  made  toward 
maintaining  laminar  flow,  and  hence  low  drag,  over  aerodynamic/ 
hydrodynamic  bodies  at  practical  flow  speeds.  In  the  case 
of  aerodynamic  bodies,  the  use  of  carefully  designed  air¬ 
foil  shapes  (e.g.  to  provide  strong  favorable  pressure  gra¬ 
dient)  and  the  use  of  boundary-layer-control  techniques 
(e.g.  suction)  has  successfully  delayed  transition  to  chord- 
length  Reynolds  numbers,  Reg^,  of  the  order  of  10  million. 

Even  larger  values  of  have  been  achieved  for  small  hydro- 

dynamic  bodies  through  the  use  of  surface  heating  and  pres¬ 
sure  gradient.  For  hydrodynamic  bodies,  extrapolations  based 
on  linear  stability  theory  indicate  that,  with  practicable 
amounts  of  surface  heating,  values  of  Rog^  in  excess  of  200 
million  may  be  possible  on  relatively  large  hydrodynamic 
bodies^ . 

The  reduction  in  drag  which  can  be  achieved  by  maintaining 
laminar  flow  over  any  vehicle  is  attractive  because  of  the 
reduced  power  requirements  to  move  the  vehicle.  However,  a 
penalty  is  generally  paid  in  maintaining  laminar  flow  on 
lifting  bodies  in  that  a  laminar  boundary  layer  separates 
much  more  easily  than  does  a  turbulent  boundary  layer,  re¬ 
sulting  in  a  significant  reduction  in  lift.  Hence,  main¬ 
taining  laminar  flow  is  practical  mainly  for  nonlifting 
bodies.  Submarines  and  torpedos  or,  more  generally,  hydro- 
dynamic  bodies  fall  into  the  latter  class. 

This  study  focuses  on  the  observed  pronounced  effects  of 
surface  heat  transfer  on  boundary  layer  transition.  Most 
Importantly,  this  report  Includes  transition  predictions 
based  on  a  relatively  new  transition  theory  for  small  and 
large  heated  hydrodynamic  bodies. 
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Section  2  presents  the  transition  equations,  including  a 
modification  needed  to  improve  transition-prediction  ac¬ 
curacy  when  surface  heat  transfer  is  present.  Included  in 
Section  3  are  transition  computations  for  various  aero¬ 
dynamic  and  hydrodynamic  boundary  layers.  The  concluding 
section  summarizes  results  and  conclusions. 


2.  FORMULATION 


2.1  EQUATIONS  OF  MOTION 

The  Saffman-Wilcox  transition  model^  ^  is  the  basic  tool 
used  in  this  study  to  analyze  effects  of  surface  heat 
transfer  on  boundary-layer  transition.  The  model's  accuracy 
has  previously  been  demonstrated  for  a  wide  variety  of  flows 
ranging  from  incompressible  boundary  layers  to  hypersonic 
blunt-body  flows.  These  applications  have  tested  the  model's 
ability  to  predict  transition  sensitivity  to  effects  of  free- 
stream  turbulence,  suction,  surface  roughness  and  pressure 
gradient.  In  all  cases,  accuracy  acceptable  for  most  en¬ 
gineering  purposes  has  been  obtained. 

For  incompressible  boundary  layers  (l.e..  for  very  small  Mach 
numbers),  the  equations  of  motion  which  constitute  the  tran- 

sition  model  are: 

Mass  Conservation 


X(nJu)  t  ^(r^v) 


(1) 


Mnmpnturn  Conservation 


9u  ^  3u  _ 


^  .  1  fpj  (u  +  -)  —1 

=  -37  ^  J  W  n' 


(2) 


Energy  Conservation 
3(C  T)  3(C  T) 


3 


Tiirhiilent  Energy 


3e  , 


pv 


(4) 


Tiirhtilent  Dissipation  Rate 


DU 


[a||fl-6p!l]pn“ 


(5) 


ir  Equations  1-5,  x  and  y  are  orthogonal  coordinates  parallel 
to  and  normal  to  a  body  surface.  The  quantity  r  is  e  ra 
coordinate  from  the  body's  symmetry  axis  while  ,1=0  for  wo 
dimensional  flow  and  J=1  for  axlsymmetrlc  flow.  The  veloc 
components  in  the  x  and  y  directions  are  denoted  by  u  an  v. 
fluid  density,  temperature,  pressure,  specific  heat,  an 
viscosity  are  denoted  by  p,  T,  p,  Cp,  and  u  respec 
quantities  e  and  (1  are  specific  turbulent  energy  and 
turbulent  dissipation  rate;  their  ratio,  e/D,  is  e  ur 
eddy  viscosity. 

Both  air  and  water  are  considered  in  this  study. 
thermodynamic  properties  for  air  are  related 
feet  gas  law  and  the  Sutherland  viscosity  law.  ppen 
lists  the  pertinent  thermodynamic  properties  of  water. 

several  empirical  parameters  appear  in  Equations 
transitional  flows,  past  studies^''  have  established  the 

following  values: 


Pr^  •■=  0.89 


e  =  0.15 
0  =  0 . 50  > 


,  =  0.09 

a*  =  0.50 


(6) 
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(7) 


a*  =  0.30  [1  -  (1  -  X*)(l  -  Re^/R^)  (1  - 

a=  .1638  [1  -  (1  -X)(l  -  Re^/R^)H(1  -  Re^/R^)]  (8) 

where  H(x)  is  the  Heaviside  r.tepfunction  and  Re^  is  turbulent 
Reynolds  number  defined  by 

Re^  =  e/fiy  (9) 

j. 

Finally,  the  parameters  \*  and  R^  are  constants  whose  values 
are 


X*  =  0.105  ,  ’’q  =  0.10 


(10) 


The  quantity  A  depends  upon  the  freestream  turbulence  level, 
T' ,  and  pressure  gradient  parameter.  A,  defined  by 


T< 


and 


Pw 


(11) 


(12) 


The  functional  dependence  of  A  upon  A  and  T'  is 

A=  .  105(1  +  2H(A)  j  1  -  exp  [-i^0A]  [  exp  [-3T'M)  (13) 

Equations  1-3  are  the  time-averaged  conservation  equations 
v'ith  classical  eddy-viscosity  and  eddy-heat-dif fusivity  clo¬ 
sure  approximations.  Equations  4  and  5  are  -onlinear  diffusion 
equations  for  e  and  fi  which  provide  a  description  of  the  growth 
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of  disturbances  in  a  laminar  flow  up  to,  through,  and  beyond 

2 

transition.  Wilcox  presents  a  thorough  discussion  of  the  way 

in  which  the  model  equations  are  used  to  predict  boundary- 

layer  transition.  Boundary  conditions  suitable  for  flow  over 

both  smooth  and  rough  surfaces  are  given  by  Wilcox  and 
h 

Chambers  . 


2.2  NEUTRAL  STABILITY  CONSIDERATIONS 


Qualitative  features  of  model-pi edlcted  transition  can  be 
conveniently  determined  by  dropping  the  convection  and  dif¬ 
fusion  terms  in  the  turbulent-energy  equation  (Equation  ^). 
In  doing  this  we  can  determine  the  neutral-stability  point 
which  is  defined  as  the  point  in  a  boundary  layer  where  tur¬ 
bulent  energy  generation,  a^^pe  |  3u/9y  |  ,  just  balances  turbu¬ 
lent  energy  dissipation,  S*p^J2e.  Hence,  neutral  stability 
is  defined  by  the  following  condition: 


=  1 

pfi  ^ 


(1^) 


Then,  noting  for  laminar  boundary  layers  that  and  a*  are 
approximately 


20  p 

6  py^ 


(15) 


a*=  0.30A^» 


Equation  1^  simplifies  to 


m^x 


y ^  |3U|  _  40 

V  Ipl  "  IT 


do) 


(17) 
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For  zero-pressure-gradient  boundary  layers  in  air,  the  velocity 
distribution  can  be  written  as 


u  =  UgF(n)  (18) 


Using  Equations  l8  and  19,  Equation  17  simplifies  to  the 
following: 

(20) 

Now  for  Incompressible  flow,  F(ri)  is  only  weakly  dependent 

upon  surface  heat  transfer.  Hence,  Equation  20  implies  that 

the  neutral  stability  Reyi  olds  number,  varies  Inversely 

as  Experiments^  Indicate  that  transition  Reynolds  number, 

2 

U  x,/v  ,  is  inversely  proportional  to  (T  /T  )  .  Assuming  that 
the  transition  Reynolds  number  is  proportional  to  the  neutral 
stability  Reynolds  number  (as  it  often  is)  implies  that  X* 
should  be  proportional  to  T^/T^. 

Computations  with  X*  Independent  of  verified  that  such  a 

dependence  is  needed.  The  unmodified  model  predicted  that 
heating  (cooling)  has  a  stabilizing  (destabilizing)  effect  on 
aerodynamic  boundary  layers,  in  contrast  to  the  measured  de¬ 
stabilizing  (.stabilizing)  effect.  These  predictions  were  un- 
surprising  as  Shamroth  and  McDonald"^  find  a  similar  reversal 
in  predicted  transition  sensitivity  to  heating  and  cooling 
with  their  turbulence-model  transition  method.  Shamroth  and 
McDonald  resolve  the  problem  by  making  a  parameter  similar  to 


X*  an  increasing  function  of  T  /T  .  While  this  model  re- 

W  0 

vision  is  sufficient  to  yield  accurate  transition  location 

for  aerodynamic  boundary  layers,  jt  is  inadequate  for  the 

hydrodynamic  case.  That  is,  according  to  Equation  20,  tran- 

2 

sition  location,  x.  ,  is  proportional  to  P,,/X*  .  For  a 

bn  W 

llqujd,  y  varies  as  T  "  near  room  temperature  (see  Appen- 

dix  /> .  ..ence,  assuminp-  X*  is  pronorrlonal  to  T..  would  ir.rlv 

rhat  X  varies  as  T  Therefore,  increasing  surface  tem- 

z  w 

perature  would  decrease  x^ ,  i.e.,  heating  would  destabilize 
a  hydrodynamic  boundary  layer.  Since  heating  stabilizes 
such  a  boundary  layer,  a  different  modification  is  clearly 
needed . 

The  parameter  X*  should  more  appropriately  depend  upon 
On  the  one  hand,  y,,/y  increases  as  T  /T  increases  for  air, 
while,  on  the  other  hand,  y,Vb^  decreases  as  T  /T  increases 
for  water.  Thus,  if  X*  were  proportional  to  some  power  of 
the  model  could  accurately  predict  effects  of  sur¬ 
face  heat  transfer  on  both  aerodynamic  and  hydrodynamic 
boundary  layer  transition.  Since  y'^1'  Tor  air,  the  fol¬ 
lowing  revised  form  of  Equation  -0  is  proposed: 

X<f  =  0.105  (b^b,,)^^^  ClOa) 

An  additional,  less  obvious,  modification  is  needed.  Note 
that  the  parameter  X  partially  control^  the  rate  at  which 
disturbances  are  amplified  beyond  the  neutral  stability 
point.  Furthermore,  if  X  ever  becomes  sufficiently  large 
(relative  to  X*)  so  that  dissipation-rate  production, 
apn^l9u/9y|,  overtakes  dissipation-rate  dissipation, 

(see  Equation  5)  before  the  neutral  stability  point  is 
reached,  transition  may  never  occur.  Thus,  if  X  remains 
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unaltered,  the  possibility  exists  that  a  finite  amount  of 
cooling  (heating)  will  cause  an  aerodynamic  (hydrodynamic) 
boundary  layer  to  forever  remain  laminar,  a  physically  un¬ 
realistic  prediction.  Hence,  Equation  13  must  be  repxaced 
by: 

X=. 105(y^/Pg)^^^  I  1  +  2H(A)[1  -  exp(-40A) ]exp(-3T*  m}  (13a) 

All  of  the  computations  presented  in  the  next  section  have 
been  performed  using  the  transition  model  defined  by  Equations 
1  through  9,  10a,  11,  12,  and  13a;  as  noted  earlier,  appro¬ 
priate  surface  boundary  conditions  are  given  by  Wilcox  and 

Chambers^ . 
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3.  APPLICATIONS 


In  the  first  part  of  this  section,  the  transition  model  is 
used  to  predict  transition  sensitivity  to  sur.’ace  heat  trans¬ 
fer  for  smooth-  and  rough-surface  aerodynamic  boundary  layers; 
computed  transition-point  locations  are  compared  with  cor¬ 
responding  experimental  data.  Then,  effects  of  surface 
heating  on  a  small  hydrodynamic  body  are  computed;  qualita¬ 
tive  comparisons  are  made  with  experimental  data.  Finally, 
transition  location  is  predicted  on  a  large  hydrodynamic 
body . 

2^1  HEATED  AND  COOLED  AERODYNAMIC  BOUNDARY  LAYERS 

One  of  the  easiest  of  all  flows  to  analyze  is  the  incompres¬ 
sible  flat-plate  boundary  layer  (FPBL).  Furthermore,  analyz¬ 
ing  this  flow  provides  a  good  test  of  the  transition  model 
as  detailed  measurements  have  been  made  to  determine  trans’'- 
tion  sensitivity  to  surface  heat  transfer.  Using  an  incompres¬ 
sible  version  of  DCW  Industries'  EDDYBL  computer  code  ,  tran¬ 
sition  computations  were  performed  for  an  incompressible  PPBL 

with 

0.5  <  <  3.0  C21) 

As  in  all  computations  in  this  study,  the  value  of  at  the 
boundary-layer  edge,  ,  was  given  by 

=  .0l85U^/yg  <^22) 

a  value  generally  used  in  EDDYBL  transition  calculations.  Re¬ 
sults  of  the  computations  are  shown  in  Figure  1;  experimental 
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data  of  Zysina  -  Molozhen  and  Kuznetsova  are  included  in 
the  figure  for  comparison. 

Because  the  freestream  turbulence  level  in  the  experiments  was 
not  available,  computations  were  done  for  a  relatively  high 
level  (T'=1.25^)  and  for  a  relatively  low  level  (T'  =  .035?).  As 

shown  in  the  figure,  freestream  turbulence  level  has  only  a 
slight  effect  on  transition  Reynolds  number  based  on  surface 
conditions,  (^^xt ^w^^e^t^’^w *  Consistent  with  the  measure¬ 
ments,  computed  is  approximately  Inversely  propor¬ 

tional  to  (T^/Tg)^;  computed  values  of  (Re^^^)^  generally  are 
within  about  20^  of  corresponding  measured  values. 

A  second,  more  subtle,  test  of  the  theory  is  to  apply  the 
equations  to  rough-wall  aerodynamic  ooundary  layers.  Surface 

cooling  can  actually  reverse  its  stabilizing  role  to  one  of 
destabilization  when  the  surface  is  rough.  The  explanation 
for  this  phenomenon  is  well  understood.  The  destabll  zation 
occurs  because  cooling  thins  the  boundary  layer,  thus  making 
the  roughness  look  larger  relative  to  a  boundary-layer  thick¬ 
ness  such  as  momentum  thickness,  6.  Hence,  the  surface  looks 
rougher,  whereby  transition  occurs  earlier. 

To  test  the  model's  ability  to  predict  this  phenomenon,  com¬ 
putations  were  performed  for  Mach  5  airflow  past  a  roughened 
spherical  body  with  nose  radius  of  2.5  inches;  roughness 
heights,  k,  between  1.5  mils  and  10  mils  were  used  in  the 
computations.  Freestream  unit  Reynolds  number  ranged  from 
3  to  10  million  per  foot.  (The  compressible  version  of  the 

2j 

model  equations  was  used  for  these  calculations).  The  ratio 

of  surface  temperature  to  freestream  total  temperature,  T.  , 

u 
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was  varied  from  0.2  to  0.8.  Figure  2  compares  results  of 
the  computations  with  a  correlation  of  experimental  data^. 
Consistent  with  the  measured  effect,  momentum  thickness 
Reynolds  number  at  transition.  Re.,  increases  with  de- 
creasing  surface  temperature. 


5  10  20 


Figure  2 .  Comparison  between  computed  and  measured 
effects  of  surface  cooling  on  rough-wall 
aerodynamic  boundary-layer  transition. 
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3.2 


HEATED  HYDRODYNAMIC  BOUNDARY  LAYERS  ON  SMALL  BODIES 


The  second  round  of  applications  is  to  a  small  heated  hydro- 
dynamic  body.  The  body,  shown  in  Figure  3,  is  a  Parsons- 
Goodson-Goldschmied®  (PPG)  minimum-drag  hull  shape  which  will 
be  referred  to  as  the  H-2  body;  the  PPG  body  parameters  are 
given  in  Appendix  B.  The  objective  of  the  computations  was 
to  determine  incipient  transition  conditions  for  various 
surface  temperatures.  By  definition,  incipienc  transition 
occurs  when  the  transition  point  is  located  at  the  maximum 
body  radius. 


All  computations  were  performed  witli  the  following  conditions 
specified: 


Ambient  Temperature,  T^ 
Roughness  Height,  k 
Freestream  Turbulence, 
Freestream  Dissipation  Rate, 


=  71°F 
=  16  pin 
=  .01f» 

=  .0185  U^/Pg 


The  pressure  distribution  was  obtained  from  the  Douglas- 
Neumann^  potential  flow  program;  Figure  shows  the  computed 
pressure  coefficient,  Cp,  defined  by 


p  =  (P-P„)/(1/2PU^) 


(23) 


In  the  figure,  z  is  axial  distance  from  the  stagnation  point. 
Three  series  of  computations  were  performed  in  which  surface 
temperature  was  held  constant  at  values  of  71  1>  76  F,  and 
8l°F.  Incipient  transition  Reynolds  number,  Re^  (based  on 
freestream  flow  conditions  and  where  V  is  body  volume), 

was  computed  for  each  value  of  T^.  Results  of  the  calculations 
are  summarized  in  Table  1;  in  addition  to  Re^,  the  table  gives 
freestream  velocity,  U„,  transition  Reynolds  number  based  on 
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Figure  4.  Computed  pressure  coefficient  for 
the  small  hydrodynamic  body 


In  Equation  24,  s  =4  feet  (the  maximum  radius  point),  and 

U 

q  (s)  Is  the  local  heat  transfer.  Figure  5  shows  computed  Re„ 
w  ^ 

as  a  function  of  T  -T  and  q^.,.  As  can  be  seen  from  the  figure, 

we  w 

heating  stabilizes  the  boundary  layer  and  hence  Increases  the 
value  of  Re^.  This  prediction  Is  qualitatively  consistent 
with  the  observed  stabilizing  effect  of  surface  heating  when 
the  fluid  Is  a  liquid.  The  magnitude  of  the  effect  Is  also 
realistic^ . 


Table  1.  Computed  Incipient  Transition  Conditions  for 
the  H-2  Body  for  Varying  Surface  Temperature. 


-  ^kwatts^ 
QwC-— r-) 


U  (knots) 

00 

10"^Re^ 

25.3 

9.4 

35.0 

13.0 

47.4 

17.6 

10"^Re< 


18.2 

24.5 

33.8 


An  Interesting  feature  of  the  predictions  Is  revealed  In 

Figure  5.  Specifically,  Re^  Increases  approximately  linearly 

with  for  the  range  of  temperatures  considered.  However, 

Re,.  Increases  much  less  rapidly  with  q  .  Hence,  Increasing 
V  ”2 

the  heating  rate  beyond  about  2  kwatts/ft  may  not  yield 

substantial  Increases  In  Re^. 
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3.3  HEATED  HYDRODYNAMIC  BOUNDARY  LAYERS  ON  LARGE  BODIES 

The  final  application  of  the  transition  model  is  to  a  large 
heated  hydrodynamic  body.  Again,  a  PPG  hull  shape,  referred 
to  as  the  R-9a  body,  is  used  (see  Appendix  B).  The  objective 
of  these  computations  was  to  determine  transition  point  lo¬ 
cation  for  specified  flow  conditions. 

Two  computations  were  performed;  in  one,  the  freestream  flow 
speed,  U^,  was  35  knots  and  in  the  other  was  30  knots.  In 
both  computations,  the  following  conditions  were  imposed: 

Ambient  Temperature,  T^ 

Surface  Temperature,  T^ 

Roughness  Height,  k 
Freestream  Turbulence,  T' 

Freestream  Dissipation  Rate, 

Figure  6  shows  the  Douglas-Neumann-computed  pressure  distribu 
tlon. 


=  55°F 

=  85°P 
=  l6  yin 
=  .01% 

=  .0185  Ug/Ug 


Tablu  2  and  Figure  7  summarize  results  of  the  computations. 
For  both  flow  speeds,  translton  occurs  well  upstream  of  the 
maximum  radius.  Computed  arclength  transition  Reynolds  num¬ 
bers,  Res^,  are  not  much  larger  than  the  value  achieved  on 
the  H-2  body  with  =  1.30  kwatts/ft^  (see  Table  1).  Noting 
the  heating  rates  Involved,  the  predictions  are  consistent 
with  results  of  the  preceding  subsection. 
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4.  SUMMARY  AND  CONCLUSIONS 


Results  reported  in  Section  3  demonstrate  the  transition 
model's  accuracy  for  boundary  layers  v;ith  surface  heat 
transfer,  provided  nonuniform  viscosity  effects  are  taken 
into  account  (Subsection  2.2).  With  no  parameter  adjustment, 
the  model  accurately  predicts  transition  Reynolds  number 
for  (a)  incompressible  smooth-wall  aerodynamic  boundary 
layers  with  surface  heating  and  cooling,  (b)  cooled  rough- 
wall  boundary  layers  on  blunt  bodies  in  a  hypersonic  air- 
stream,  and  (c)  heated  hydrodynamic  boundary  layers  on  small 

bodies . 

Results  of  the  hydrodynamic  computations  agree  with  the 
measured  strong  boundary-layer  stabilization  attending 
small  amounts  of  surface  heating.  The  model  also  predicts 
that  large  amounts  of  heating  are  not  significantly  more 
effective  in  delaying  transition  than  the  smaller  rates 
considered  in  the  H-2  body  computations  of  Subsection  3.2. 
The  larger  R-9a  body  computations  support  this  prediction. 
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APPENDIX  A 

THERMODYNAMIC  PROPERTIES  OF  WATER 

The  following  thermodynamic  properties  of  water  are  perti¬ 
nent  to  boundary  layer  transition:  mass  density  (p),  speci¬ 
fic  heat  (C  ),  thermal  conductivity,  (k),  molecular  viscosity 
(y),  \nd  laminar  Prandtl  number  (Pr^).  Values  used  in  this 
study  for  these  quantities  have  been  obtained  from  Schlicht- 
ing^^and  are  valid  for  temperatures  ranging  from  about  hO°F . 
to  110°F. 

Mass  density,  specific  heat,  and  thermal  conductivity  are 
approximately  constant  and  have  the  following  values: 


=  1.936 

Ibf  •sec^/ft‘* 

(Al) 

=  1 

Btu/lbf ‘“F 

(A2) 

=  0.35 

Btu/ft ‘hr*  °P 

(A3) 

Viscosity  and  Prandtl  number,  by  contrast,  are  strongly  tem- 
peratui'e  dependent  over  this  range  of  temperatures.  An  ap¬ 
proximate  polynomial  fit  to  the  Schlichting  data  was  used  to 
calculate  y  in  the  computations  of  Section  3,  l.e., 

-6  —8 

y  =  7.943*10  (T/600)  Ibfsec/ft^  (a4) 

with  temperature,  T,  in  degrees  Rankine .  Finally,  the  laminar 
Prandtl  number  follows  from  its  definition,  namely,  PrL=yCp/K, 
wherefore , 

Pr.  =  2.51(T/600)"®  (A5) 

jj 

Figure  A1  compares  Equations  A4  and  A5  with  the  measured  de¬ 
pendences  of  y  and  Prj^  upon  temperature. 
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APPENDIX  B 

HYDRODYNAMIC  BODY  SHAPE  PARAMETERS 


The  two  hydrodynamic  bodies  analyzed  in  this  study  are 
specified  in  terms  of  the  eight  dimensionless  parameters 
used  by  Parsons,  Goodson,  and  f’oldschmied  .  The  values  of 
these  parameters  are  as  follows: 

Parameter  H-2  Body  R-9a  Body 


f 

4.450 

5.000 

r 

""m 

0.427 

0.600 

k. 

1 

2.715 

i^.OOO 

0.260 

0.050 

^i 

0.740 

0.700 

1.650 

1.750 

0.595 

0.800 

t 

0.180 

0.150 

Additionally,  the  H-2  body  is  9-457  feet  long  and  has  a  volume 
of  14.5  cubic  feet;  the  R-9a  body  is  55  feet  long  and  has  a 
volume  of  2500  cubic  feet. 
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